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It  ABSTRACT 

Using  spatial  coherency  estimates  of  Rayleigh  and  Love  waves 
recorded  at  ALPA,  LASA  and  N0RSAR,  it  is  demonstrated  that  for  a 
particular  frequency  the  wave  propagation  cannot  be  described  by 
a  unique  wavcvcctor.  Instead,  the  propagation  is  better  described 
by  a  distribution  in  the  wavenumber  domain  or,  in  physical  terms, 
a  distribution  in  azimuth  and  phase  velocity.  In  general  the 
conerency  in  the  mean  direction  of  propagation,  which  is  a  function 
of  the  velocity  distribution,  stays  high  relative  to  the  coherency 
along  the  mean  wavefront,  which  is  a  function  of  the  azimuthal 
distribution. 
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INTRODUCTION 


Large  arrays  of  seismometers  allow  us  to  study  how 
seismic  signals  vary  over  an  area  of  the  earth's  surface. 
Some  idea  of  the  nature  of  this  variation,  and  its 
magnitude,  is  necessary  for  array  design  and  subsequent 
data  processing.  One  measure  of  this  variation  is  the 
coherence  of  the  time  shifted  waveform  between  seismo¬ 
meter  pairs.  This  paper  presents  experimental  results 
for  the  variation  of  long  period  surface  waves  at  the 
three  large  arrays  ALPA,  LASA  and  NORSAR. 

Previous  seismic  coherence  studies  have  been 
confined  to  long-period  noise  (Capon,  1969),  short- 
period  isotropic  noise  (Aki,  1957,  and  Backus,  et  al, 
1964),  and  short-period  directional  noise  (Bungum, 
et  al,  1971).  However,  the  earth  is  not  homogeneous  and 
so  the  waveform  of  a  'pure*  signal  changes  as  it  propa¬ 
gates,  due  to  refraction,  diffraction,  dispersion,  and 
scattering.  This  change  of  waveform  causes  a  loss  of 
coherence  which  depends  on  the  seismometer  separation 
and  the  wave  frequency.  There  is  a  further  coherence 
loss  due  to  the  additive  background  noise  but  in  this 
study  only  signals  with  large  s ignal-to-noise  ratios 
are  used  and  so  the  cohe  ence  loss  caused  by  the  back¬ 
ground  is  considered  to  be  negligible. 


SPATIAL  COHERENCE  AND  PROPAGATION 


When  a  surface  wave  propagates  in  a  medium  which 
has  lateral  inhomogeneities  comparable  in  scale  to  the 
wavelength,  the  constant  phase  surfaces  of  the  wave  are 
neither  simply  planar  nor  cylindrical.  The  wave  fronts 
cannot  be  represented  by  a  single  wave  vector  at  each 
frequency.  For  a  time  window  of  particular  length  used 
in  analysis,  the  wave  vector  at  each  frequency  has  a 
distribution  both  in  magnitude  and  direction,  or,  in 
other  words,  a  particular  frequency  component  of  the 
wave  does  not  have  a  unique  phase  velocity  or  arrival 
azimuth.  The  wave  is  then  represented  in  wavenumber 
space  as  a  distribution  F(k)  rather  than  a  singularity 
<S(li"k0)  at  a  particular  frequency.  The  spatial  coherence 
y(f,r)  is  then  defined  as: 


Y(f.r) 


1 

FT 


F(k,f  )exp(-2irik»r)dk 
“'K 


The  integration  is  taken  over  the  area  K,  which  is  the 
region  in  the  wavenumber  domain  where  F(k)  i  0.  N  is 
a  suitable  normalizing  factor. 

Consider  a  signal  recorded  at  two  sites  separated 
by  a  distance  r.  The  coherence  between  the  two  at  a 
frequency  f  is  normally  defined  as: 

y3  (f) .  mil 

12  G1(l)C2(t) 
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where 


G 1 2 ( ^ )  =  cross  power  spectrum  at  frequency  f 

G1 Cf )  =  power  spectrum  at  site  1 

G 2 ( f 3  *  power  spectrum  at  site  2 

By  estimating  the  coherence  of  the  signals  between 
all  seismometer  pairs  in  an  array,  the  spatial  variation 
of  the  coherence  y(f,r)  can  be  obtained,  where  r  is  the 
vector  separation  of  any  two  positions  within  the  array. 
The  wavenumber  distribution  F(k)  responsible  for  the 
spatial  coherence  can  then  be  obtained  by  inversion  of 
expression  (1).  If  the  distribution  F(k)  is  confined  to 
a  relatively  small  region  of  wavenumber  space,  it  has 
been  demonstrated  theoretically  (Gossard,  1969)  and 
experimentally  (Mack  and  Flinn,  1971)  that  the  two- 
dimensional  transformation  of  expression  (1)  can  be 
well  approximated  by  independent  one-dimensional  trans¬ 
formations,  with  the  result  that  the  variation  in 
coherence  in  the  direction  of  propagation  is  a  function 
of  the  wavenumber  magnitude  range  (and  hence  the 
velocity  range)  whereas  the  variation  in  coherence 
normal  to  the  direction  of  propagation  is  only  a  function 
of  the  angular  distribution  of  F(k).  To  illustrate  this 
consider  the  wavenumber  geometry  shown  in  Figure  1. 

The  wavefield  is  represented  by  an  area  of  finite 
thickness  2Ak  subtending  an  angle  ±  0  about  the  origin. 
This  wavenumber  geometry  describes  a  wave  group  propa¬ 
gating,  on  the  average,  in  the  negative  Y  direction.  The 


words  "on  the  average"  are  used  because  the  azimuth  of 
propagation  is  really  distributed  ±  0  about  the  Y  axis. 
Similarly  the  magnitude  of  the  wavevector  is  distributed 
as  kQ  ±  Ak.  Physically,  this  means  that  the  phase 
velocity  has  a  distribution  rather  than  a  unique  value 
at  a  particular  frequency,  within  the  area  of  the  array. 
The  range  of  velocity  is  given  by: 


(v,,  v,) 


r  f 

[ftat  * 

o 
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If  only  small  spread  angles  and  velocity  ranges  are 
considered,  as  was  mentioned  before,  the  coherence  loss 
in  the  wavefront  direction  is  virtually  independent  of 
the  loss  in  the  wave  propagation  direction.  This  is 
equivalent  to  making  the  region  K  rectangular,  in  which 
case  the  integral  (1)  becomes  the  product  of  two  simple 
integ  als. 


kQsin0 


■  o 

=  7k  sln5  F(k  ,f)exp(-2  i 

°  J-k  dnfl 


kQsin0 


ikxx  dkx 


1 

•  m 


k0>Ak 


k0-Ak 


F(ky,f)exp(-2rikyy)dky 


For  various  simple  forms  of  F(k) ,  equation  (4)  can 
easily  be  solved.  The  four  forms 
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1. 


F(kx)  •  1,  F(ky)  -  1 

2.  F(kx)  •  exp(-a|kx|),  F(ky)  -  exp(-B|kv  -  kQ|) 

3.  F(k/  «  exp(-ak2) ,  F(ky)  -  exp(-B(ky-k0) 2 

4.  F(k)  •  fiCk-kj)  ♦  6(k-k2) 

imply  spatial  coherence  of  the  forms: 

sin(2irk  xsinS)  sin(2*Aky) 

1.  y(r,f)  ■  - 2 -  .  - 

2Tikoxsin0  2irAky 


2.  Y(r,f) 


2a  28 

“2 - 7  *7 - 7 

a  +4ttx  8  *4uy* 


3.  Y  (r ,  f) 


4 •  Y(r,f)  -  [1+cos 2tt (IC2-1CJ ) »jr ]/2 


These  expressions  have  the  property  that  for  y»0 
the  coherence  parallel  to  the  wavefront  is  a  function 
of  the  angular  scatter,  and  s.milarly  for  x-0  the 
coherence  normal  to  the  wavefront  is  a  function  of  the 
velocity  scatter. 

Model  (1)  is  perhaps  physically  unreasonable  in  the 
sense  that  it  is  difficult  to  imagine  a  scattering  pro¬ 
cess  where  the  scattered  waves  have  the  same  amplitude 
within  a  certain  angle  and  zero  everywhere  else.  Some 
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sort  of  decay  with  increasing  angle  is  more  to  be 
expected.  Models  (2)  and  (3)  satisfy  this  condition. 
However,  there  is  not  much  difference  numerically  in 
the  three  models  for  given  ranges  of  scattering  and  all 
three  have  similar  shapes,  i.e.,  an  increasingly  negative 
first  derivative  in  the  region  of  small  separation. 

Models  (2)  and  (3)  are  suggested  by  Chernov  (1960) 
for  wave  propagation  in  random  media.  Model  (4)  is  the 
expression  which  is  representative  of  the  multipathing 
phenomenon  described  by  Capon  (1969).  Two  or  more 
discrete  waves  associated  with  the  same  phase  (e.g.,  a 
Rayleigh  wave)  propagate  across  the  array  and  so  exhibit 
two  distinct  wavevectors.  This  last  form  causes  the 
most  severe  variations  in  spatial  coherence. 

Figures  2  and  3  show  theoretical  curves  for  coher- 
2 

cnce  (y  )  versus  seismometer  separation  for  the  case  of 
Model  (1)  at  a  period  of  21.3  seconds.  In  Figure  2  the 
separation  is  taken  along  the  mear.  wavefront,  and  the 
various  curves  illustrate  how  the  coherence  varies  with 
increasing  angular  scatter.  Figure  3  shows  coherence 
versus  separation  in  the  direction  of  propagation  for 
various  velocity  ranges.  Figure  4  and  5  show  the  spatial 
coherence  behavior  at  a  period  of  25.6  sec,  again  for 
Model  (1).  These  particular  periods  are  used  in  the 
theoretical  examples  because  with  the  window  length  used 
in  the  spectral  analysis  these  two  discrete  frequencies 
lay  in  the  part  of  the  spectrum  having  the  most  power. 
Figure  6  illustrates  the  effect  of  multipathing,  using 
Model  (4),  on  spatial  coherence.  Using  such  standard 
curves  for  comparison  with  coherence  measured  in  the 
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two  orthogonal  directions  allows  the  principal  dimensions 
of  the  wavenumber  distribution  to  be  estimated  and  hence 
the  angular  scatter  and  velocity  scatter  of  the  wave. 


Effect  of  dispersion 

In  estimating  the  coherency  of  a  wave  train  between 
two  points  it  might  be  assumed  that  dispersion  would 
contribute  to  the  loss  of  coherence  in  the  direction  of 
propagation.  The  phase  relationship  within  a  band  of 
frequencies  would  vary  with  time,  hence  propagation 
distance. 

This  would  cause  changes  in  the  real  and  imaginary 
parts  of  the  Fourier  transform  which  are  smoothed 
separately  in  the  coherence  estimation.  Using  the  phase 
velocity  dispersion  curve  for  LASA  given  by  Glover  and 
Alexander  (1969),  a  synthetic  example  was  used  to 
determine  the  dispersion  contribution  to  the  coherency 
loss.  The  result  was  that  for  a  distance  separation  of 
100  km  the  coherence  was  still  as  high  as  0.99  for  a 
period  of  20  seconds.  It  will  be  seen  that  the  observed 
coherence  loss  is  significantly  greater  than  this  so 
it  is  reasonable  to  assume  that  dispersion  contributes 
very  little  to  the  form  of  the  observed  spatial  coherency 
at  LASA.  It  has  been  assumed  that  the  result  would  be 
similar  at  NORSAR  and  ALPA. 


COHERENCE  ESTIMATES  OF  RAYLEIGH  WAVES 

Estimates  have  been  made  of  the  spatial  variation 
of  the  coherence  of  surface  waves  at  LASA,  NORSAR  and 
ALPA.  The  signals  were  aligned  before  the  transformations 
and  subsequent  spectral  smoothing  were  performed.  Align¬ 
ment  is  necessary  in  order  to  eliminate  the  phase  shift 
caused  by  the  relative  delays  across  the  array.  Smoothing 
a  spectrum  which  includes  such  a  phase  shift  leads  to 
an  erroneous  value  for  the  cross  power  spectrum  and  hence 
an  incorrect  coherency  estimate.  Twenty-four  degrees  of 
freedom  were  used  in  the  coherency  estimation  resulting 
in  acceptably  narrow  confidence  intervals  at  the  90% 
level . 

NORSAR 

Two  examples  using  NORSAR  recordings  are  presented 
to  illustrate  the  azimuthal  dependence  of  spatial 
coherence.  Figure  7  shows  the  spatial  coherence  of  a 
Rayleigh  wave  from  an  earthquake  in  North  Si.nkiang,  the 
individual  channels  being  shown  in  Figure  8.  The  period 
under  consideration  is  21.3  seconds}  and,  although  there 
appears  to  be  a  small  separation  in  the  two  orthogonal 
directions,  the  overall  coherence  remains  high  across 
the  maximum  extremities  of  the  array.  Comparison  of  the 
measured  results  with  standard  curves  reveals  that  an 
angular  scatter  of  ±  4°  about  the  mean  azimuth  explains 
the  loss  of  coherence  along  the  wavefront.  The  slight 
decrease  of  coherence  in  the  direction  of  propagation 
can  be  explained  in  terms  of  a  phase  velocity  scatter 
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of  approximately  ♦  0.1  km/sec  about  the  mean  value. 

The  spatial  coherence  of  the  Love  wave  at  the  same 
period  is  shown  in  Figure  9.  Again  the  coherence  remains 
high  across  the  array  and  has  about  the  same  distribution 
as  the  Rayleigh  wave. 

At  other  azimuths  the  picture  can  be  decidedly 
different.  Figures  10  and  11  show  a  Turkish  event 
recorded  at  NORSAR  and  the  spatial  coherence  of  the 
Rayleigh  wave  respectively.  The  signal-to-noise  ratio 
is  high  and  the  signal  appears  visually  to  have  little 
variation  across  the  array.  However,  the  spatial  coher¬ 
ence  falls  off  much  more  rapidly  with  sensor  separation 
than  in  the  previous  case.  The  difference  in  coherence 
in  the  two  orthogonal  directions  is  now  very  obvious. 

In  this  situation  lines  of  equal  coherence  are  roughly 
elliptical  with  the  major  axis  pointing  in  the  mean 
direction  of  propagation.  If  the  angular  distribution 
is  continuous,  as  in  Models  (1-3),  a  scatter  of  +_  10° 
about  the  mean  direction  explains  the  loss  of  coherence 
along  the  wavefront.  If  the  variation  is  attributable 
to  two  discrete,  equal-amplitude,  interfering  waves, 
represented  by  Model  (4),  then  the  azimuthal  separation 
is  about  12°. 

The  loss  of  coherence  in  the  direction  of  propaga¬ 
tion  is  somewhat  scattered  but  a  standard  curve  for  a 
velocity  range  of  +0.3  km/sec  about  a  mean  value  lies 
through  the  measured  values.  If  this  coherency  loss  is 
caused  by  two  discrete  waves,  or  modes,  propagating 
from  the  same  back  azimuth  the  phase  velocity  difference 
is  approximately  0.25  km/sec. 


-9- 


ALPA 

The  Sinkiang  event  used  in  the  NORSAR  analysis  was 
also  analyzed  using  the  ALPA  recordings.  The  individual 
channels  and  beam  are  shown  in  Figure  12  and  the  spatial 
coherence  of  the  Rayleigh  wave  at  a  period  of  21.3  sec 
is  shown  in  Figure  13.  It  is  obvious  that  the  loss  of 
coherence  along  the  mean  wavefront  is  quite  severe  and 
corresponds  to  a  continuous  angular  distribution  of  t  20° 
or  to  two  discrete  waves  separated  by  30°. 

Frequency  wavenumber  analysis  of  this  particular 
Rayleigh  wave  showed  the  primary  arrival  crossing  the 
array  with  a  back  azimuth  of  314°.  A  second  Rayleigh 
wave,  presumably  multipathed,  crossed  the  array  with  a 
back  azimuth  of  about  350°  almost  at  the  same  time. 

This  would  certainly  explain  the  rapid  decrease  in 
coherence. 

The  scatter  in  the  coherence,  even  though  the 
signal-to-noise  ratio  is  very  high,  suggests  that  the 
scattering  is  more  complex  than  can  be  explained  by  just 
one  of  the  models. 

Presumably  a  'discrete'  multipathed  arrival  has  a 
finite  distribution  in  the  wavenumber  plane  so  the 
measured  spatial  coherence  would  be  affected  by  both. 

The  loss  of  coherence  in  the  direction  of  propaga¬ 
tion  is  rather  scattered  but  appears  to  decrease  in  a 
consistent  saanner  with  increased  separation.  A  curve 
passing  through  the  points  indicates  that  a  velocity 
scatter  of  about  ±  0.35  km/sec  about  a  mean  value 
could  explain  this.  However  this  value  may  be  too  high 
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because  the  assumption  that  the  coherence  is  independent 
in  the  two  orthogonal  directions  starts  to  be  invalid 
for  wavenumber  distributions  which  subtend  angles 
greater  than  20°  (Gossaid  and  Sailors,  1970). 

In  an  attempt  to  estimate  how  much  scattering  occurs 
in  the  Alaskan  region  itself,  the  Rayleigh  waves  generated 
by  Cannikin  on  Amchitka  Island  and  recorded  at  ALPA  were 
subjected  to  the  coherency  analysis.  The  results  are 
shown  in  Figure  14.  The  separation  :n  the  two  orthogonal 
directions  is  again  quite  striking.  Even  though  the 
epicentral  distance  is  only  about  20°  the  loss  of 
coherence  along  the  wavefront  suggests  scattering  over 
a  range  ±  12°  about  the  mean  azimuth.  The  coherence 
remains  very  high  in  the  direction  of  propagation  which 
indicates  that  the  phase  velocity  at  this  period  is 
almost  single  valued  or,  in  other  words,  there  is  very 
little  velocity  scatter  or  mode  mixing. 

The  lower  spatial  coherency  values  at  ALPA  with 
respect  to  NORSAR  for  Central  Asian  events  helps  to 
explain  the  constant  surface  wave  magnitude  difference 
for  the  same  event  measured  at  both  arrays  (Mack,  1972). 
The  marked  difference  in  coherence  in  the  two  directions 
suggests  that  elliptical  arrays  with  the  major  axis 
pointing  towards  the  source  would  give  the  best  signal- 
to-noise  ratio  improvement.  This  has  been  demonstrated 
for  some  Rayleigh  waves  recorded  a+  ALPA. 

LASA 

The  two  events  used  to  investigate  spatial  coherence 
at  the  LASA  also  highlight  the  difference  between  the 
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scattering  of  a  'pure'  Rayleigh  wave  and  a  multipathed 
example.  Figure  IS  illustrates  the  latter  situation.  In 
Figure  15a  the  orthogonal  estimates  01  coherence  are 
shown  for  a  window  400  seconds  long  and  in  Figure  15b 
this  has  been  extended  to  700  seconds.  Frequency  wave- 
number  analysis  showed  the  initial  Rayleigh  wave 
arrival  from  a  back  azimuth  of  300°  -  305°  and  a 
second  arrival  about  four  minutes  later  from  a  back 
azimuth  of  315°  -  320°.  This  latter  wave  was  contained 
in  the  long  window  but  not  in  the  short  window  during 
the  coherency  estimation.  The  coherence  loss  for  the 
short  window  can  be  explained  by  an  azimuthal  distri¬ 
bution  of  about  ♦  5°  and  a  velocity  of  ±  0.1  -  0.2  km/sec. 
about  the  mean  value. 

The  rapid  fall  off  in  coherence  using  the  longer 
window  suggests  a  second  discrete  arrival  separated  by 
about  20°  from  the  first  and  this  is  in  agreement  with 
the  frequency-wavenumber  analysis.  It  can  be  seen  that 
the  second  arrival  has  had  no  effect  on  the  velocity 
dependent  coherence. 

The  second  LASA  example  is  for  a  Rayleigh  wave 
from  an  event  in  Baja  California.  The  propagation  path 
is  entirely  continental.  Inspection  of  the  individual 
channels  in  Figure  16  indicates  good  signal  similarity 
and  no  apparent  interference.  However  the  spatial 
coherence  does  show  separation  in  the  two  orthogonal 
directions  and  this  can  be  explained  by  an  azimuthal 
range  of  15°  about  the  mean  and  a  velocity  range  of 
t  0.2  km/sec. 
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DISCUSSION 


Spatial  coherence  estimates  at  array  sites  help  us 
to  understand  the  propagation  of  surface  waves  in  a 
heterogeneous  medium. 

In  addition  to  the  phenomenon  of  raultipathing 
described  by  Capon  (1970)  it  would  appear  that  ’pure' 
surface  wave  phases  do  not  exhibit  a  single  propagation 
direction  or  velocity.  The  unique  wave  vector  has  to  be 
replaced  with  a  distribution,  the  dimensions  of  which 
can  be  estimated  approximately  from  the  coherence.  The 
wavenumber  distribution  estimated  in  this  way  is  inde¬ 
pendent  of  the  array  response.  However,  if  the  diameter 
of  the  array  is  too  small,  the  functional  form  of  the 
spatial  coherence  is  left  in  doubt  because  all  the 
models  suggested  in  this  paper  have  approximately  the 
same  shape  for  the  shorter  separation  distances  and 
there  is  sufficient  scatter  in  the  estimates  so  that 
all  the  theoretical  curves,  or  combinations  of  them, 
tend  to  fit  the  measured  estimates  of  the  coherence.  In 
this  light,  it  is  obvious  that  the  quantitative  values 
given  for  the  wavenumber  distribution,  reflecting  the 
scattering,  must  be  considered  as  order-of-magnitude 
estimates. 

The  azimuthal  variation  strongly  indicates  that  at 
least  part  of  the  scattering  is  a  function  of  the  travel 
path  rather  than  the  structure  under  the  array. 

Tne  wave  scattering  phenomenon  decreases  with 

increase  of  period  and  is  very  small  for  periods  greater 
than  40  seconds. 
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The  fact  that  the  phase  velocity  at  periods  of 
about  20  seconds  has  a  distribution  rather  than  a 
single  value  limits  the  resolution  with  which  this 
portion  of  the  dispersion  curve  can  be  inverted  to 
obtain  crustal  structure.  This  distribution  is  in 
the  range  0.1  -  0.4  km/sec  even  for  'pure*  Rayleigh 
waves . 
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Figure  1.  Wavenumber  representation  of  a  wave  group  propa¬ 
gating  with  a  range  of  azimuth  and  velocity. 
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Figure  2.  Theoretical  coherence  versus  seismometer  separa¬ 
tion  along  the  mean  wavefront  and  parametric  in  range  of 
azimuth  of  the  propagation  vector  for  a  period  of  21.3  sec 
(velocity  3  4.0  km/sec). 
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Figure  3.  Theoretical  coherence  versus  seismom«-.er  separa¬ 
tion  in  the  mean  direction  of  propagation  and  parametric  i 
the  range  of  phase  velocity  for  a  period  of  21.3  sec. 


Figure  4.  Theoretical  coherence  versus  seismometer  separa¬ 
tion  along  the  mean  wavefront  and  parametric  in  range  of 
azimuth  of  the  propagation  vector  for  a  period  of  25.6  sec 
(velocity  *  4.0  km/sec). 


Figure  5.  Theoretical  coherence  versus  seismometer  separa 
tion  in  the  mean  direction  of  propagation  and  parametric 
in  the  range  of  phase  velocity  for  a  period  of  25.6  sec. 


Figure  7.  Coherency  at  a  period  of  21.3  sec  versus  seismometer 
separation  for  the  Rayleigh  wave  shown  in  Figure  8. 
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Figure  9.  Coherency  versus  seismometer  separation  for  the  Love 
wave  lrom  North  Sinkiang  recorded  at  NORSAR. 


SEISMOMETER  SEPARATION  (h«) 


Figure  10.  Coherency  at  a  period  of  21.3  sec  versus  seismometers 
separation  for  the  Rayleigh  waves  shown  in  Figure 
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Figure  15a.  Spatial  coherency  estimates  on  the  first  400  sec  of 
Rayleigh  wave  from  an  earthquake  in  the  Andreanoff  Islands.  Date 
30  March  1971;  origin  time,  11:30:38.9,  A  =  43°. 


7.  Coherency  at  21.3  sec  period  versus  seismometer  separation 
Rayleigh  waves  shown  in  Figure  16. 


